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Abstract 

The two-gluon exchange model of the pomeron is used to compute the pho- 
toproduction reaction 7p Z^p. The predicted cross section is too smah 
to be observed at HERA, but may be detectable at an eventual Next Linear 
Collider. 



I. INTRODUCTION 



The two-gluon exchange model [0-0] has for a long time offered a semi-quantitative 
understanding of the pomeron in QCD. The model originated as a view of elastic scattering. 
It has been applied extensively to vector meson photoproduction, and has recently been 
used to predict jet production in double pomeron interactions [0,0]. In this paper, we use 
the model to calculate exclusive Z° photoproduction at low Q^. 

The diagrams for 7p Z^p are shown in Figure 1. The upper half of each diagram can be 
calculated in a straight-forward manner, since only known electroweak and gluonic couplings 
appear. The lower half of each diagram involves the non-perturbative color structure of the 
proton, but the required discontinuity can be obtained by parametrizing it and constraining 
the parameters to fit pp elastic scattering. The calculation of photoproduction thus 
introduces no new parameters. A similar calculation can be done for 7p — > Tp, where 
provides a large momentum scale; but that calculation is not as clean because it brings in 
the T hadronic wave function. 

Reggeization of the gluons and interactions between them (e.g., as given by ladder di- 



agrams) must be responsible for the gradual energy dependence of pomeron exchange [jlO 



We will ignore these effects here. From a practical standpoint, the energy dependence of 
— > Z^p will be dominated by the minimum longitudinal momentum transfer that is 
kinematically required when the energy is not far above threshold. 

The predicted cross section is necessarily small, because the large mass mz sets a short 
distance scale ~ l/rriz for the typical transverse separation between the qq pair that form 
the Z^. This leads to a factor l/m'^z ^^e final amplitude, which is proportional to the 
color dipole moment of the pair. A couple of enhancement factors oc hrniz arise in the 
calculation, as explained below; but the final predicted cross section turns out to be too 
small to observe at the HERA ep collider in its present form. 
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II. CALCULATION OF 7P ^ Z^P 



The imaginary part of the amphtude is equal to l/2i times the discontinuity given by 
the four diagrams of Figure 1. The upper half of diagram (a), for example, corresponds to 
a factor 

J (27r)4 [{p, - - m2] [{p, - k^f - m2] " ^ ^ 

in the discontinuity, where 

tr„ = rr{7-e(pi) [m-i'^-{ki - Pi)] 7/. ("i + ^7-^2) (A + ^75) 

'^■e*{p^)[m-i'^-{ki- q2)]'^y{m-i'-fki)} . (2) 

To obtain the high energy limit of the amplitude, we use light-cone coordinates p± = {po ± 
Pz)/V2 where pi+ and p2- are large and s = pi^p2- . One of the two delta functions in 
Eq. ([1|) is saved for the eventual integration over gluon momentum gi, which carries no large 
+ or — component as a result of this and a similar delta function from the lower half of 
the diagram. The other delta function reduces Eq. (|T]) to a three-dimensional integral over 
transverse momentum and light-cone momentum fraction x = 

Vj"^ = 5(pi-gi + 0(l))pi,pi^T„ 

T 1 r 2, P2f^P2u tig , . 

''^"47rV Jo x{l-x) l'^"'^ d^d2 ' ^ ^ 

where 

di = (pi — ki)'^ — = [{ki — xpi)\ -|- — x(l — 

d2 = (P3 - k2f -m^ = [{ki - g2 - a^ps)! + - x{l - x)pl]/{l - x) . (4) 

For = we can choose a gauge in which e{pi) has only transverse components. Then 
e(p3) can be assumed to have only transverse components as well, because the production 
of longitudinal Z° can be shown to be suppressed by a factor \/~i/mz in the amplitude. 
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Upon summing over all four diagrams, the contribution from 75 terms is found to vanish 
in the large s limit. Combining the denominators in Eq. using the Feynman parameter 
trick ^ = Jq dy [ay + 6(1 — y)]^'^ allows the ki± integral to be performed, and leads to 

Vap'-'^ = 6{pi-qi + 0(l))pi„pi^T(i) (5) 
T« = F«(g2±,A^)-F«(0,A^) (6) 

where 

F^'\q±,A^) = - f'dx ['dy{e^ip,)-elips)i[m'-yil-yy]/c 

TT JO Jo 

— [1 — 2x{l — x)] Inc) + 4:V-e±{pi) f ■£^(^3) x {1 — x)y {1 — y)/c} 
c = m'^ + y{l-y)v'^ -x{l-x) [ypj + (1 - y)pi] 

V = q± + xA± . (7) 

Here pi = —Q^ is approximately zero, p| = m|, and A_l = (^3 — pi)^ = (gi — q2)± is the 
transverse momentum transfer. 

The F{q2±, A_i_) term in Eq. @ comes from the two "off-diagonal" diagrams in which 
one gluon is attached to the quark line and the other to the anti-quark. The -F(0, A_|_) term 
comes from the diagonal diagrams in which both gluons are attached to the same line. There 
is a strong cancellation between these two contributions when qij_ or q2± is small, because 
the two-gluon system is color neutral and the qq system is spatially compact due to the large 
Z° mass. However, the diagonal term will dominate in the full amplitude because the range 
of eventual integration over is controlled by mz- 

We can estimate the integral in Eq. (|^) using the approximation m| ^ qfj_, g2±) A^^. 
This limit is obtained by splitting the integration into six separate regions according to 
< y < yo or yo < y < 1 and < x < xq, xq < a; < 1 — xq, or 1 — xq < a; < 1. The result is 
independent of Xq and yo in the limit xq -C 1, yo ^ 1 and is given by 
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+ [gi±-e±(Pi)g2±-el(p3) + g2±-e±(pi)gi±-el(p3)]/m2} 



/(a) = a [ln(m /m^) — m h &ln 
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1 + 4/a (8) 

We want the limit of high energy at small momentum transfer, where the helicity non-flip 
amplitude dominates. We can therefore take e±{pz) = e_L(pi) = ^(=i=l, —2,0,0) to obtain 

e±(pi)-el(p3) = l, 

gi±-e±(pi) g2±-el(p3) + g2±-e±(pi) gi± -61(^3) = gi±-g2± ■ (9) 

It is a good approximation to neglect the quark mass m, since Eq. (H) is not singular in the 
limit m^O. This approximation is reasonably good numerically even for quark masses up to 
the charm quark mass, in the important regions of q\^, g^x ^'^d A^, because the corrections 
to it do not contain the large factor \n.m\. In this way, we find 

g{q^) = q^[\nq^/ml-nx]. (10) 

Eq. ([T0| ) was derived for the region where and = (gi — ^2)! are small compared 

to m|. We can provide it with approximately correct behavior when these quantities become 
comparable to m| by inserting an additional factor [1 + + glx)/""^!]"^ i'^^o Eq. (|T0|). 
The adequacy of the resulting approximation to T^^^ , for the purpose of calculating the full 
amplitude, has been checked by numerical integration. 
The coefficient A in Eq. ( ]T0| ) is given by 

4 TT^ a 20 . o ^ , , , 

A = —n ^ (1 - IT sm^ 0^) ■ 11 

sm cos u,,, 9 



This includes a factor 2 from the sum over quark flavors u + d and s + c, and a factor 1/2 
from the sum over quark colors. I take = 0.25 for the strong coupling at the low average 
momentum transfer scale that occurs here. This is slightly optimistic, since the amplitude 
receives contributions from a spectrum of gluon transverse momenta extending all the way 
up to ~ niz- It leads to A = 0.083 . 

A reasonable parametrization for the discontinuity of the gluon-proton amplitude repre- 
sented by the bottom half of Figure 1 has been given in Ref. [§]: 

V^, = S{p2-qi + 0(1)) P2mP2.T(2^ (12) 

where 



with a = ^Jq± + ^m?q and h = ^(A — q)]_ + Am^ . A selection of reasonable parameter 
choices that have been tuned to fit pp elastic scattering is given in Table I. Dependence of 
the results on the choice of parameter set provides an estimate of systematic errors. 
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TABLES 
Table I 

Parameters and predicted cross section a and slope B for 7p — > Z^p. 
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N 




5[GeV-2] 


0.14 


0.5 


5.797 


2.65 X 10^ 


0.014 


7.0 


0.14 


0.3 


4.944 


4.89 X 10^ 


0.017 


6.4 


0.30 


0.5 


6.659 


2.13 X 10^ 


0.018 


7.5 


0.30 


0.3 


5.901 


1.10 X 10 "- 


0.020 


6.9 


1.0 


0.5 


7.638 


4.06 X 10*^ 


0.045 


7.9 


1.0 


0.3 


7.098 


9.04 X 10^ 


0.047 


7.5 
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The full amplitude for 7J9 Z^p in the two gluon exchange picture of Figure 1 (including 
a factor 8 from the sum over gluon colors) is given by 

with = — (pi — = —t. A finite gluon mass irig is included in the propagators to 
suppress contributions from long distance, as an approximation to color confinement [|l],0|. 
In elastic scattering, this is necessary to avoid unphysical behavior of the elastic slope in the 
limit t — > 0. The integral can be done numerically at each momentum transfer Aj^. The 
result is enhanced by two powers of lnm|: one coming directly from T*^^) (see Eq. ([lOD ) and 
the other coming from the integration over g^, since the integrand behaves as dq\/q]_ for 
large q\ up to ~ Tn^z). 

The differential cross section is given by 

da 1 



dt IGtts^ 



\M\' (15) 



in the high energy limit. It is independent of s in that limit. The predicted integrated 
cross section a = J ^ dt is listed in Table I for each choice of parameters. The result is 
approximately 0.025 pb, with an uncertainty of about a factor of 2 coming from the range 
of parameters that provide a plausible description of elastic scattering. (It is interesting to 
note that this calculated result is not far from the crude estimate of ~ 0.08 pb that can be 
made by scaling the observed cross section ioi '^p ^ J /ijjp aX HERA [|lT| by the dipole size 
factor {mji.^/mzY.) 

The dependence of ^ on momentum transfer t is roughly exponential, so we can charac- 
terize it by a slope parameter B such that ^ oc e^*. The value of B, defined by reproducing 
a and ^ at t = 0, is also listed in Table I. It is approximately equal to the one half the slope 
of elastic scattering, since it arises essentially from the wave function effect associated with 
the proton remaining intact, while in pp elastic scattering, both protons must remain intact. 
Equivalently, one may say that the extent of the interaction in impact parameter space is 
set by that of the proton. 



III. PREDICTION FOR EP EZ^P 

One can hope to look for photoproduction in electroproduction ep eZ^p, which 
will be dominated by low transverse photons producing Z°'s with the same helicity as 
the photon. The cross section is given by 

"^'^Zq^^''^ = My^^') -rviW) , (16) 

where 

'l + {l-yf 2mly 



(17) 



y Q' 

is the flux of transverse photons |[Tl|]. Defining the four- momenta in ep ^ e Z^ p as k, p, 
k', q', p' respectively, the usual kinematic variables are s = {k + p)^, = —{k' — kY, and 
y = {k — k') ■ p/k ■ p . The role of s in Sect. |I| is now played by W"^ = {k — k' + pY = ys. 
The minimum value of is given by Q^in ~ ^ly'^/i^ ~ y)- The maximum occurs 
when the electron is scattered through 180°, but most of the cross section comes from the 
region of very small Q^. (We have assumed the limit — ^ in calculating the cross section 
for — » Z^p in Sect. |T[ but if that condition is relaxed, the model shows no dramatic 
dependence on Q"^.) 

We assume that the dependence on is dominated by the effect of the minimum mo- 
mentum transfer |tmin|- Specifically, we assume the approximate dependence ~ c^* with 
B ^ 7GeV~^, which was found in the large W limit, to be approximately valid at finite W 
where the minimum value of t = {p — p'Y is different from because of the longitudinal 
momentum transfer necessary to produce the high-mass state. This assumption might ap- 
pear optimistic because producing the on-shell intermediate states in Figure 1 involves an 
increase in mass of the states at both ends of each gluon line, which requires a noticeably 
larger minimum longitudinal momentum transfer. However, the scale of transverse momen- 
tum transfer for each gluon exchange extends all the way to 0{mz), so the longitudinal 
momentum transfer associated with the individual gluon exchanges need not be dominant. 
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Integrating Eq. (|T6D at the HERA energy y/s = 300 GeV leads to (Jep^ezOp/ o-jp~*zOp = 
0.032. Combining this with the results listed in Table I yields a predicted cross section for 
diffractive production at HERA of approximately 1 femtobarn. A realistic experimental 
search would have to rely on one of the clean decay modes —>■ e+e^ or Z^ —>■ fi~^fj^~ , which 
would suppress the rate by a branching fraction 0.067. The resulting predicted rate is too 
small by 2 - 3 orders of magnitude to be observed at the HERA ep collider in its present 
form. 



IV. CONCLUSION 

The quasi-elastic photoproduction process •yp — > Z^p has been calculated using the 
two-gluon exchange model of the pomeron. The calculation is rather clean because the 
electroweak part of it is completely well defined and the soft proton part can be normalized 
to elastic scattering. 

The process would lead to Z^ production in ep collisions, which would have a clean 
and dramatic signature: one e'^e'^ or fj^^fi~ pair with invariant mass equal to the Z^ mass, 
and hence with very large individual transverse momenta; and nothing else in the detector, 
since the scattered electron and proton would generally disappear down the beam pipes. 
Unfortunately, the predicted rate is too small to be observed with the luminosity currently 
available at HERA. Possible upgrades to HERA that would increase its luminosity by 2 -3 
orders of magnitude could make the process observable. 

Because of the higher available energy, the related diffractive process 77 Z^p'^ should 



be observable with a 77 option at a (9(1 TeV) Next Linear Collider [jT2[. The cross section 
can be estimated from Sect. || and Regge factorization as ~ 0.3 femtobarns for 77 energies 
far above threshold. 
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